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ABSTRACT 

We investigate the energetics of the jet and X-ray corona of Cyg X-l. We show that 
the current estimates of the jet power obtained from Ha and [O in] measurements 
of the optical nebula surrounding the X-ray source allow one to constrain the bulk 
velocity of the jet. It is definitely relativistic (v >0.1c) and most probably in the range 
(0.3-0.8)c. The exact value of the velocity depends on the accretion efficiency. These 
constraints are obtained independently of, and are consistent with, previous estimates 
of the jet bulk velocity based on radio measurements. We then show that the X-ray 
emission does not originate in the jet. Indeed, the energy budget does not allow the 
corona to be ejected to infinity at relativistic speed. Rather, either a small fraction of 
the corona escapes to infinity, or the ejection velocity of the corona is vanishingly low. 
Although the corona could constitute the jet launching region, it cannot be identified 
with the jet itself. We discuss the consequences for various X-ray emission models. 

Key words: accretion, accretion discs - black hole physics - ISM: jets and outflows 
- radiation mechanisms: non-thermal - X-rays: binaries - radio continuum: stars 



1 INTRODUCTION 

Cyg X-l is the prototype of black hole binaries. It is a per- 
. sistent X-ray source, powered by accretion onto a black hole 
from a massive companion, HDE 226868, most likely via a fo- 
cused wind. The value of the mass of the black hole is subject 
to controversy, it is in the range of M ~ (5-15)Mq accord- 
ing to Herrero et al. (1995) or M ~ (14-27)M (Ziolkowski 
2005; see also Gies & Bolton 1986). In this paper we will 
assume M — 15Mq. The distance, D, is most likely within 
D ~ 2.1 ± 0.2 kpc (Ziolkowski 2005 and references therein); 
hereafter we adopt D — 2 kpc. 

As other black hole binaries Cyg X-l presents various 
X-ray spectral states. The main ones are the so-called Low 
Hard State (hereafter, LHS) and the High Soft State (here- 
after HSS, see Done, Gierlihski & Kubota 2007 for a com- 
prehensive review). In the LHS, the isotropic bolometric lu- 
minosity Z/h,obs is rather stable and close to 2xl0 37 erg s -1 
(see e.g. Gierlinski et al. 1997; Frontera et al. 2001; Di Salvo 
et al. 2001, McConnell et al. 2002; Cadolle Bel et al. 2006). 

Cyg X-l is also a variable radio source. Multi- 
wavelength observations in the LHS have shown the presence 
of a flat-spectrum radio emission which is produced by syn- 
chrotron emission from relativistic electrons in a compact, 
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self-absorbed jet (Hjellming & Johnston 1988). This jet was 
later resolved at the milliarsec scale by Stirling et al. (2001). 
The average radio flux is about 15 mjy in the LHS. In the 
HSS this radio emission is strongly suppressed, indicating 
that the radio jets may not be produced. Deep radio obser- 
vations of the field of Cyg X-l resulted in the discovery of a 
shell-like structure which is aligned with the resolved radio 
jet (Gallo et al. 2005). This large-scale (5 pc in diameter) 
structure appears to be inflated by the inner radio jet. Gallo 
et al. (2005) estimate that in order to sustain the observed 
emission of the shell, the jet of Cyg X-l has to carry a kinetic 
power that is comparable to the bolometric X-ray luminos- 
ity -Lh.obs of the binary system. Then Russell et al. (2007) 
refined this estimate using Ha and [O ill] measurements of 
the jet-powered nebula. They estimate that the total kinetic 
power of the double sided jet is Lj = (0.9-3) xlO 37 erg s — 1 . 
If we adopt Lh=2 x 10 37 erg s" 1 as the typical X-ray lu- 
minosity in the hard state then j = Lj/Lh is in the range 
0.45-1.5. 

The bolometric luminosity does not change dramati- 
cally during the state transition to the soft state (see e.g. 
Frontera et al. 2001; Zdziarski et al. 2002; Malzac et al. 
2006; Wilms et al. 2006). In the following we combine this 
fact with the estimate of the jet power to constrain the ini- 
tial and terminal jet velocities, as well as the nature of the 
X-ray emitting region. The structure of the paper is as fol- 
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Figure 1. Constraints on the jet terminal velocity as a function 
of the measured jet power. The conservative lower limit on /3oo as 
a function of j is given by the black thick solid curve (assuming 
the maximum possible value for A = 4, and minimum value for the 
accretion efficiency in the soft state r] s = 0.06 and an efficiency of 
1 in the hard state) . The other curves show the dependence of fScxj 
on j in the more reasonable case in which the efficiency remains 
constant across the state transition. The red hatched region shows 
the range of /3oo allowed for A = 3 and the accretion efficiency r/ s 
ranging between 0.06 and 0.4. The blue hatched region shows the 
same for A = 2 . The full lines show the dependence of /3 oo on j 
for the fiducial case r] B = 0.1. The black thin solid lines show the 
upper and lower limits provided by the radio imaging and radio 
X-ray correlations (Gleissner et al. 2004). The black thin dotted 
lines show the constraints from the modelling of the super-orbital 
periodicity (Ibragimov et al. 2007). The horizontal grey stripe 
shows the overlapping region between those two constraints. 



lows. In section we write down the equations for the en- 
ergy budget of the accretion/ejection flow and show that the 
accretion efficiency cannot increase much during the state 
transitions from LHS to HSS. Then in section |2~21 we show 
that under reasonable assumptions the terminal jet velocity 
must be relativistic. Then in section [2731 we show that if the 
terminal jet velocity is indeed relativistic and the Thomson 
depth along the radius of the base of the jet is larger than 
unity, then the initial velocity of the plasma in the base of 
the jet must be non- relativistic (and probably very low). 
These results and their caveats are then discussed in the 
context of current accretion models in section [4] 



2 ENERGETICS OF CYG X-l 

2.1 Energy and mass budget, accretion efficiencies 

We define Mh as the mass accretion rate in the hard state 
and at sufficiently large distances from the black hole so that 
it is representative of the amount of material available for 
both accretion and ejection. A fraction fj of the accreting 
mass is ejected, the rest is swallowed by the black hole. The 
total power output of the system (jet+radiation) must be 
equal to that extracted through accretion: 



where rjh is the accretion efficiency in the hard state, and c 
the speed of light. The accretion efficiency characterises the 
total energy available for both radiation and jet production. 
The radiative efficiency can be expressed as : 

L b 1 - fj 



M h c' 1 



= Vh- 



1+3 



(2) 



Therefore if we knew both the radiative and total accretion 
efficiency, the known estimate of the jet power would al- 
low us to evaluate the ejected mass fraction /j. In turn, the 
ejected mass fraction would yield the jet velocity. However 
both efficiencies are poorly constrained, they depend on the 
unknown and possibly complex dynamics of the accretion 
flow in the hard state (possible role of advection, possible 
non-keplerian orbits). The form of the accretion flow in the 
hard state is still a matter of debate. In order not to rely 
on any specific model we will consider that 7/h may take any 
value comprised between and 1. 

Then, in order to obtain some meaningful constraints 
we will have to combine this with the simpler situation that 
occurs in the HSS. Indeed, in the soft state there is obser- 
vationally no evidence for a jet. We will therefore assume 
that the jet is not produced in the HSS, or at least that it is 
energetically negligible. This means that, in the soft state, 
the radiative efficiency and the total accretion efficiency are 
identical: 



■ T] s M s c 



(3) 



where L s , r; s and M s are respectively the luminosity, accre- 
tion efficiency and mass accretion rate in the HSS. There is 
overwhelming evidence that accretion in the soft state pro- 
ceeds predominantly through a geometrically thin disc. The 
accretion efficiency of such a disc is theoretically limited in 
the range r) a =0.06-0.4 depending on the black hole spin. 
Then, using, equations ([2]) and ([3]), we define: 



A 



(4) 



L] + — 7jh(l — fj)MhC 



(1) 



L B M h _ 77s 1 + j 
L h M s % 1 - /j ' 

A is nothing else than the ratio of the HSS to LHS radiative 
efficiencies. The value of this ratio is a crucial parameter de- 
termining the energetics of the system. It is reasonably well 
constrained by the observations. In the HSS, the isotropic 
bolometric luminosity is somewhat higher than in the LHS. 
Observations performed at different epochs with different 
instruments provide estimates of the isotropic bolometric 
luminosity in the HSS, L s ,obs ^ (6.2-7.2) xlO 37 erg s -1 
(Gierlinski et al. 1999; Frontera et al. 2001; McConnell et al. 
2002). Therefore, during the state transition the observed 
luminosity jumps at most by a factor Z/ s ,obs/£h,obs — 3- 
4 (Zdziarski et al. 2002). We note that owing to possible 
anisotropy of the X-ray emission, the observed isotropic lu- 
minosities may be different from the real intrinsic luminosi- 
ties entering the definition of A. However, as will be shown 
in section |31 the effects of radiation beaming are expected to 
be weak and usually tend to reduce A. Moreover the tran- 
sition to soft state is likely to be triggered by an increase 
in mass accretion rate and therefore Mh/M s < 1. For these 
reasons, we can safely constrain A ^ 4. 

If, as it is widely believed, accretion proceeds in an ad- 
vection dominated accretion flow in the hard state (Narayan 
& Yi 1994; Esin et al. 1997) then one would expect to have 
77s > 77h. In fact, equation Q implies 77s/77h < A/(l + j). As 



© 2009 RAS, MNRAS 000, [T]{9] 



Energetics of Cyg X-l 3 



we know that A ^ 4 and j ^ 0.45, this gives rjs/rfb ^ 2.75. 
Therefore, the accretion efficiency does not increase dramat- 
ically across the transition from LHS to HSS. Moreover, 
since ??h must be lower than unity, and r/ s is larger than 
0.06, we can set very conservative limits on the efficiency 
ratio: 0.06 sC Tj s /r] h sC 2.75. 

In the following we will use A = 3, j = 1, r) s = rfa = Q.l 
as typical values. ?7 S =0.1 corresponds to a black hole with a 
moderate spin parameter a = 0.4. We note however that the 
recent spectroscopic results of Miller et al. (2009) suggest a 
nearly non spinning black hole (a < 0.06) in which case the 
efficiency r/ a would be close to 0.06. We will also consider the 
unlikely combination A = 4, r/ s — 0.06, r/^ = 1, j=0.45 as an 
extreme set of parameters used to provide a robust limit on 
the jet velocities. 

Finally, we can write the fraction of ejected material as: 

i+i 



fj = i 



It. 



A 



(•») 



For the fiducial model the fraction of ejected material is 
fj — 1/3. For the extreme parameters most of the accreting 
material is ejected (fj — 0.978). 



2.2 Terminal jet velocity 

Based on the absence of detection of the counter jet, Stirling 
et al. (2001) give a lower limit on the bulk velocity of the 
radio jet of poo > 0.6. Similar considerations and the lack 
of response of the radio emission on short time-scales led 
Gleissner et al. (2004) to constrain the jet velocity in the 
range 0.4 < Poo < 0.7. Ibragimov, Zdziarski and Poutanen 
(2007) find a similar result 0.3 < Poo < 0.5 from the mod- 
elling of the super-orbital modulation observed in the X-ray 
and radio band. In the following we compare these estimates 
with the one provided by the jet and disc energetics. 
The jet kinetic power is: 



Lj = jL h = f 3 -M h (joo - l)c 



(6) 



where is the 'terminal' Lorentz factor of the jet at large 
distance from the black hole, in the region where most of 
its power is deposited in the interstellar medium. Combin- 
ing equations and ([5]), it can expressed in terms of the 
accretion efficiencies: 



7=0 = 1 + 



A -(.7 + 1); 



(7) 



For a typical accretion efficiency r/ s = rjh = 0.1, assum- 
ing j = 1 and A = 3, we derive Poo ~ 0.4, in remarkable 
agreement with the constraints based on the radio observa- 
tions. Fig. [1] explores how this result depends on the effi- 
ciency parameters. 

The red hatched area of Fig. [T] shows the region of the 
Poo-j plane that is allowed by the uncertainty on the ac- 
cretion efficiency rj a (for A = 3 and rfc = rj s ). This region 
overlaps widely with the constraints from the radio obser- 
vation and allows Poo =0.2-0.9. 

As mentioned above A could be lower than 3, if so the 
velocity of the jet must be larger. As shown on Fig. (TJ for 
A = 2, (3oo is in the range 0.3-1. For A = 2 (and rj s = Tjh) 
there is no solution with j > 1 (see equation @). 

Fig.[2]shows that the terminal jet velocity increases with 
the ratio r/ s /rjh and may become very large. Indeed for a 
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Figure 2. The jet terminal velocity as a function of of the ac- 
cretion efficiency ratio J/s/^h f° r a jet power of 3xl0 37 erg s — 1 
(top panel) and 9 xlO 36 erg s _1 (bottom panel). In both panels 
the blue, red and green curves stand for A=2,3 and 4 respectively. 
The dashed, full and dotted curves stand for ?7s=0.06, 0.1 and 0.4 
respectively. The curves are plotted only in the range of f/s/^h 
for which there is a solution (/j > 0) and 77^ < 1. The black thin 
solid lines show the upper and lower limits provided by the radio 
imaging and radio X-ray correlations (Gleissner et al. 2004). The 
black thin dotted lines show the constraints from the modelling 
of the super-orbital periodicity (Ibragimov et al. 2007). The hor- 
izontal grey stripe shows the overlapping region between those 
two constraints. 



given j and A there is a critical value of the ratio ?7s/?/'h for 
which the ejected fraction fj vanishes (see equation [5]) and 
the jet terminal velocity becomes infinite. There is no so- 
lution above this critical value of r/s/^h- As seen on Fig. [2] 
for ,7=1.5 the maximum possible rj s /rjh are 0.8, 1.2, 1.6 for 
A=2,3 and 4 respectively. For a weaker jet, the critical values 
of rja/r/h are increased. For j = 0.45, they are 1.38,2.05 and 
2.75 for A=2,3 and 4 respectively. This confirms that the ac- 
cretion efficiency cannot increase much during the transition 
from LHS to HSS. 

Overall the terminal bulk velocity increases with rj a , j, 
and rjs/rih and decreases with A. The dependence of Poo on 
j in the extreme case r/s = 0.06, 77h = 1, A = 4 is shown in 
FigE If provides a conservative lower limit on the terminal 
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Figure 3. Constraints on the jet initial velocity and radius in 
Cyg X-l. The hatched region shows the exclusion area delimited 
by the relation /3o vs r$ in the extreme case A = 4, ri s = 0.06, 
jjh = 1, j = 0.45, corresponding to a minimal terminal velocity 
/3oc = 0.12. The region above the thin continuous line is allowed in 
principle but requires the jet to slow down (/3o > /3oo)- The dashed 
line shows the relation between r and /3rj for j = 1 a terminal jet 
velocity (3^ =0.4 which is favoured by the radio observations and 
also corresponds to our fiducial model with A = 3, r} s = »7h = 0.1. 
In all curves the optical depth r = 1 and the black hole mass is 
15 M . 

jet velocity /3oc > 0.12. We conclude that, if the nebula of 
Cyg X-l is powered by the jet kinetic power, the terminal 
velocity of the jet must be at least mildly relativistic. 

2.3 Size and velocity of the X-ray emitting region 

Taking into account the fact that the jet is two-sided, mass 
is ejected at a rate: 



Mj = 2-KR 'yof3ocn mp 



(8) 



where m p is the proton mass, no is the comoving frame 
density at the base of the jet, R$ its radius, 70 and f3oc 
are the Lorentz factor and velocity of the material entering 
the jet. Then, combining equations and © and setting 



L h =2xl0 37 erg s~\ we find: 



roJoPo = 2 x 10" 



J 



15M 







M 



(9) 



T(7oo - 1) 

where r = nourRo is the Thomson depth along the radius 
of the base of the jet, ro = Rq/Rq and Rq is the gravita- 
tional radius. Assuming there is no jet in the soft state and 
that the optical nebula is powered by the jet kinetic power, 
equation (|9]) may be combined with equation (0 to express 
the initial jet velocity in terms of the accretion efficiencies 
parameters rj s and «h and the observationally constrained 
parameter A. 

For now, let us assume that the X-ray emission in the 
hard state is produced at the base of the jet. If so r should 
correspond to the Thomson depth of the X-ray emitting 
medium. Let us now further assume that the X-ray emission 
is produced by Comptonisation. This standard assumption 
is overwhelmingly supported by the data (see e.g. Zdziarski 
& Gierlihski 2004). Then, using Comptonisation models, the 



Thomson depth r can be measured directly through spectral 
fits. 

There is a continuum of spectral parameters between 
the hard and soft state rather than sharp transition. To de- 
fine the spectral states one needs to set some boundaries on 
the spectral parameters. The exact value of these boundaries 
are somewhat arbitrary. It is customary to consider that the 
source is in the hard state when the X-ray photon index 
T is in the range 1.5-2.1. But the 'softer' hard state spec- 
tra are less common and appear usually when the source is 
about to make a transition or failed transition (e.g., Malzac 
et al. 2006). These softer spectra are called 'hard interme- 
diate states' by some authors (Homan & Belloni 2005; Del 
Santo et al. 2008). They are probably not representative of 
the stable hard state spectra associated with steady compact 
radio jets. For this reason in this paper we will consider the 
parameters when the source is far away from the transition 
i.e. when the spectrum is really hard F ~ 1.6. 

When the source is clearly in the hard state, spectral 
fits with Comptonisation models yield a Thomson optical 
depth tt in the range 1-3 and electron temperature fcT e ~ 
60-100 keV (see e.g. Gierlinski et al. 1997; Frontera et al. 
2001; Cadolle Bel et al. 20060 depending on the details of 
the models and the observation. It is also apparent in the 
fig. 6 of Ibragimov et al. (2005) and fig. 5 of Wilms et al. 
(2006) that the hardest spectra of Cyg X-l are concentrated 
in a small region of the TT-kT e plane with tt ^ 1. 

Comptonisation models used to fit the data assume a 
corona with a simple geometry such as a sphere or a slab. 
The Thomson optical depth is defined along the smallest 
dimension of the corona (i.e. the radius for the sphere, or the 
height of the slab). Only the smallest dimension is relevant 
since photons escape preferentially along this direction. In 
the case of emission in the base of the jet, the geometry 
will be cylindrical. If the scale height of the X-ray emitting 
region is larger than Rq then we should have r = tt — 1— 
3. If on the contrary the vertical extension of the corona is 
smaller than R$ we simply have r >l-3. In any case, the 
optical depth r entering in equation (|9} is larger than unity. 

Then, as a consequence of the large optical depth of 
the base of the jet, either the corona is tiny (r < 4r 9 ) or 
the velocity of the corona is very small. Indeed, for a given 
r, setting f3oo = 0.12 (i.e. to the minimum value allowed by 
the considerations of section [2T2]| in equation Q provides an 
upper limit on the initial jet velocity. Fig. [3] shows that for 
t = 1 this limit implies that for any reasonable jet section 
(ro — 10 — 100) the initial velocity must be non-relativistic 
(Po <0.1). In the far more likely case of a jet terminal ve- 
locity /3oo = 0.4, that is consistent with both the previous 
radio estimates and the present ones, we find that f3o must 
be lower than a few 10 -2 for any reasonable r. 

We conclude that, if the Comptonising corona consti- 
tutes the base of the jet (or more generally if the base of 

1 An earlier spectral analysis using EXOSAT data suggested the 
optical depth could be lower tt ~ 0.3 and the temperature larger 
kT c ~ 150 keV (Haardt et al. 1993). Nevertheless, more recent, 
better quality data confirm the optical depth is large. For tt < 1 
the individual Compton scattering orders become apparent in the 
high energy spectrum, producing bumps that are not observed. 
Also the higher temperature does not provide a good fit to the 
very sharp high energy cut-off. 
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the jet has a Thomson depth r > 1) then its vertical bulk 
velocity is non-relativistic. 



3 BEAMING EFFECTS 

In our estimates of the ratio of luminosities of the two spec- 
tral states, we have so far neglected the anisotropy of the 
X-ray emission. In the HSS this anisotropy is due to a geo- 
metric effect: the emission is produced in the thin disc plane. 
In the LHS, the X-ray emitting region has a velocity /3o lead- 
ing to Doppler beaming of the radiation. This could affect 
our estimates of A and the constraint on the jet velocities. 
In order to take these effects into account we rewrite A as: 



2.0 



A 



M h 



fb, 



(10) 



where L Si0 b s and Lh.obs are the isotropic luminosities ob- 
served in the soft and hard state respectively, fb is a factor 
accounting for beaming effects and anisotropy. We know that 
the ratio ~(3-4) and that |k < 1. 

We have now to estimate /b. Since in the soft state 
the luminosity is predominantly emitted as blackbody radi- 
ation, and since L s is the power emitted by the two-sides 
of the disc, we have L s , bs — 2L s /i (with /i = cosi and i is 
the angle between the line of sight and the normal of the 
accretion disc). Then following Beloborodov (1999) for the 
Doppler beaming effects of the X-ray emitting region with 
bulk velocity /3o, fb can be evaluated as follows: 



/b 



1+X 



(11) 



4 7o V(i-/W 

Because we consider the total energy radiated by both faces 
of the disc, this expression differs from the formula given by 
Beloborodov (1999) by a factor 1/2. The factor \ accounts 
for additional sources of radiation in the hard state (other 
than the direct emission from the corona). For instance, in 
the case of an accretion disc corona above a cold thin disc 
some level of disc reprocessing and reflection is expected. 
Then, assuming isotropy at the disc surface: 



1 + /V2 



(12) 



l + /3o 2 ' 

Alternatively, in the absence of a cold disc \ may ac- 
count for the X-ray radiation coming from the opposite side 
of the accretion flow and transmitted through the corona. 
Then in the optically thin limit: 



1 + /?oM / 

On the basis of analysis of absorption lines Gies & 
Bolton (1986) estimated i=33±5 degrees. The polarimetric 
measurements of Dolan & Tapia (1989) yield i in the range 
25-67 degrees. The upper panel of Fig. [3] shows that in the 
allowed range of inclinations the fb factor is of order of unity 
except for large flo when the Doppler beaming actually de- 
creases significantly fb. Therefore the effects of anisotropy 
will tend to increase the value of /3oo rather than the oppo- 
site. For instance we can set A = 4/b, r/ a = 0.06, ?7h = 1 and 
j = 0.45 and compute the resulting minimum /3oo as a func- 
tion of /3q. The middle panel of Fig.|4]shows that at small /3o 
our lower limit is not significantly affected and stays in the 
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Figure 4. Effects of beaming and anisotropy. Top panel: depen- 
dence of the beaming factor /j, on the velocity of the X-ray emit- 
ting region /3q. Middle panel: effects on the minimum terminal 
velocity, /3oo as a function of fio- In the parameter space below 
the diagonal dotted line the terminal jet velocity is lower than 
the initial jet velocity. Bottom panel: effects of beaming on the 
maximum initial jet velocity. /3q as a function of r. In all panels 
the blue curves stand for an inclination of 25 degrees, black for 
33 degree, and red for 67 degree. The solid lines assume a corona 
above an accretion disc (making use of equation 11121 ) while the 
dashed lines dashed represent the results for a corona without 
optically thin disc and taking into account the radiation from the 
opposed side of the corona (using equation J 13t ") ■ The thick black 
curves are for the fiducial model A = 4/b, rj s — T]h = 0.1 and 
j = 1. While the other thinner curves are for the limiting case 
A = 4/b, r; s = 0.06, r] h = l and j = 0.45. 
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range 0.1-0.2. If the initial velocity approaches the speed of 
light the terminal velocity must be large. We also plot the 
Pao-Po relation expected for the more likely set of parame- 
ters rj s = 0.1, ?7h = 0.1 and j = 1, in this case the terminal 
jet velocity must be > 0.3c, again in agreement with the 
constraints from the radio data. 

As shown in section 12.31 a lower limit on the terminal 
jet velocity sets an upper limit on the initial velocity at the 
base. The bottom panel of Fig. [3] shows how the constraints 
discussed in section [2T3l mav be affected by beaming. Keeping 
A = 4/b, 77s = 0.06, r/h = 1 and j = 0.45, the upper limit on 
0o as a function of r is not changed qualitatively. Moreover 
except for the largest possible inclination angle, beaming 
tends to reduce /3o. For the most likely inclination of 33 
degree, the upper limit on /3o is reduced by a factor of ~ 2. 
For the fiducial set of parameter, rfa = rjh = 0.1 j = 1 and 
an inclination of 33 degree, j3o < 0.01 for any reasonable r. 



4 DISCUSSION 

4.1 Consequences for the jet structure 

4-1.1 Dark jet component, bending and misalignment 

Assuming that the jet is quenched in the soft state, and that 
the optical nebula around Cyg X-l is powered by the jet 
kinetic energy, we have set a lower limit on the terminal jet 
velocity in the hard state of Poo > 0.1. The actual velocity 
is most likely in the approximate range 0.3-0.8. In section|3l 
we investigated the possible effects of beaming of the X-ray 
radiation and found that they do not affect qualitatively 
these conclusions. 

Our estimates of /3oo are the first constraints on the 
jet velocity obtained independently of any the radio mea- 
surements. This is important because Gleissner et al. (2004) 
and Ibragimov et al. (2007) constrain only the velocity of 
the radio emitting material. There could be additional jet 
components which do not produce any radiation. Indeed, 
modelling the extended radio emission observed by Stirling 
et al. (2001), Heinz (2006) finds that the kinetic power of 
the radio jet is several orders of magnitude lower than the 
estimates obtained from the observations of the optical neb- 
ula. In order to explain this puzzling result, Heinz (2006) 
suggests several alternatives. Among these he discusses the 
possibility that the source of the kinetic energy powering the 
nebula is not the steady radio-emitting jet but some other 
dark component. Our results indicate that if there is a dark 
component its velocity is comparable to that of the radio 
emitting material. 

Modelling the radio orbital periodicity in terms of free- 
free absorption in the wind of the companion, Szostek & 
Zdziarski (2007) also infer the presence of a second jet 
component, which must be not only dark but also slow 
(Ax> ~ 5 x 10 _2 c). Indeed, in their model, the observed phase 
lags between radio bands can be explained if the jet is bent. 
As shown by these authors, the bending may be caused by 
a slow dark component. This slow component could trans- 
port most of the material and energy. Our results are not 
consistent with this picture and require an average velocity 
which is definitely relativistic. Perhaps, the jet is bent for a 
different reason, such as the effects of ram pressure from the 
wind of the companion star. 



We also note that even if the jet is not bent, similar 
phase lags can be produced by free-free absorption (as in the 
model of Szostek & Zdziarski 2007), provided that the jet 
average direction is not exactly perpendicular to the orbital 
plane but tilted by a few degrees. This would not be uncom- 
mon among X-ray binaries. Indeed, Maccarone (2002) re- 
ports that in at least two microquasars (GRO J 1655-40 and 
SAX J 1819-2525) the observed relativistic jets appear mis- 
aligned. This is believed to be caused by the spin of the black 
hole not being perpendicular to the orbital plane. Then the 
central parts of the accretion disc are forced to rotate in the 
same plane as the black hole (Bardeen & Petterson 1975). 
The jets are produced by the inner part of the disc, per- 
pendicular to it and therefore misaligned. But then the disc 
also exerts a torque on the black hole which will eventually 
align the black hole spin (see e.g. Natarajan & Pringle 1998; 
King et al. 2005). However the current theoretical estimates 
indicate that the alignement time-scales are likely to be at 
least a substantial fraction of the lifetime of these systems 
(Maccarone 2002). For Cyg X-l, a rough estimate can be 
obtained from equation 2.16 of Natarajan & Pringle (1998) 
which gives an alignment time of order of 4xl0 5 yr (for a 
spin parameter a =0.06). Unfortunately, the time since the 
formation of the black hole is poorly constrained. It could be 
comparable to the optical nebula lifetime, i.e. (2-6)xl0 4 yr 
according to Russell et al. (2007) , possibly much longer, but 
shorter than the time since the formation of the progenitor, 
i.e. (3.5-6.5) x 10 6 yr according to stellar evolution models 
(Mirabel & Rodrigues 2003). If the black hole was formed 
recently with a misaligned spin, the jet of Cyg X-l could 
still be misaligned nowadays. This may be the cause of the 
radio lags. 

4-1-2 Jet acceleration and mass loading 

In section [231 it was shown that if the Comptonising corona 
is being ejected to infinity (or more generally if the base of 
the jet has a Thomson depth r > 1) then its vertical bulk 
velocity is non-relativistic. In the case of a Compton corona, 
the observation requires a Thomson depth larger than unity 
along the vertical scale height of the corona. Therefore, the 
velocity must be low over a vertical scale h which is at least 
of order of r/r. Otherwise, if acceleration occurred at a lower 
scale height, the corona would be depleted of its material 
and be optically thin. Then, the jet material leaving this 
launching region has to be gradually accelerated at larger 
distances in order to reach the mildly relativistic speed ob- 
served at large scales. 

It is not clear however which acceleration mechanism 
could produce such a velocity profile. In magnetically driven 
jet models most of the acceleration occurs relatively close 
to the accretion disc (see e.g. Ferreira 1997; Casse & Fer- 
reira 2001). In principle, more distant acceleration could be 
achieved by converting the gas internal energy into bulk ki- 
netic energy through longitudinal pressure gradients. If most 
of the initial thermal energy can be converted into bulk ki- 
netic energy, and neglecting other forms of internal energy, 

the terminal velocity fiao is of order of ./ 6 ' a , where 2] is 

the gas temperature at the base of the jet. In consequence 
accelerating the jet up to mildly relativistic velocities re- 
quires a temperature that is a few times virial: /3oo > 0.3 
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would imply T\ > 2 x 10 11 K. Such high temperatures are 
expected in two temperature accretion flows (Narayan & 
Yi 1994). Similarly, Merloni & Fabian (2001) elaborate on 
the possibility of launching strong thermally driven outflows 
from 2-temperature accretion disc coronae. We note however 
that in the case of Cygnus X-l, the ion temperature of the 
hot corona is likely to be too low. Indeed, Malzac & Bel- 
mont (2009) show that the relatively large Thomson depth 
implies Ti < 2 x 10 10 K which appears too cold to drive a 
strong thermal outflow. 

Perhaps a more likely situation would be that only a 
small fraction of the X-ray emitting material is loaded onto 
the jet and ejected to infinity. This would allow a Thomson 
depth larger than unity in the corona and arbitrarily low in 
the base of the jet. The initial velocity could then be rela- 
tivistic, depending on the fraction of ejected material. This 
would also avoid the requirement for most of the jet acceler- 
ation being produced at large distances from the black hole. 
This scenario would also be consistent with recent numerical 
simulations (e.g. Hirose et al. 2004; McKinney 2006; Fragile 
& Meier 2009) suggesting that the jet launching region is 
compact, initial launching velocity is relativistic, with low 
mass loading onto the jets. 

4.2 Consequences for accretion/ejection models 

4-2.1 X-ray jet models 

When the importance of the connection between radio and 
X-ray emission was realised, it was proposed by several au- 
thors that the hard X-ray emission of the LHS could orig- 
inate from the jet. Most of these models require or assume 
an initial velocity that is too large. In order to reproduce 
the observed X-ray luminosity, the external Compton model 
of Georganopoulos, Aharonian & Kirk (2002) requires a jet 
power of ~ 5 x 10 38 ergs which is immediately excluded by 
the estimates of Russell et al. (2007) . In the model by Reig, 
Kylafis Giannios 2003 (see also Giannios, Kylafis & Psaltis 
2004; Kylafis et al. 2008) the soft photons from the accre- 
tion disc are upscattered by the jet. In this model the optical 
depth is in the range r = 1.2 — 15, and the radius of the base 
r = 75 - 300, with the velocity (3 ~ 0.5 - 0.8. This is in 
clear breach of the constraints of section 12.31 

These conclusions do not apply to the jet model of 
Markoff et al. (2001, 2005 herafter M05). In this model the 
optical depth at the base of the jet is very small (r <C 1) 
and makes it energetically possible to have a mildly rela- 
tivistic initial jet velocity. Indeed, M05 fit the spectral en- 
ergy distribution of Cyg X-l with a jet dominated model 
which attributes the X-ray emission to thermal synchrotron 
self-Compton in the base of the jet. These fits result in a jet 
power that is comparable to the X-ray luminosity and there- 
fore consistent with the estimates of Gallo et al. (2005) and 
Russell et al. (2007). M05 do not quote explicitly the result- 
ing optical depth at the base of the jet but it can be easily 
estimated from the other parameters: they infer a radius at 
the base ro which is in the range the 4-10 (depending on the 
data set). The initial jet velocity is "foPo — 0.4 and the ter- 
minal jet velocity ~ 2-3. Using equation ((9]l this implies 
an optical depth r in the range 10" 3 - 3xl0" 2 . This is much 
lower than what is usually found when the same data are 
fit with thermal comptonisation models (i.e. r > 1). In fact, 



in this model the optical depth is so small that the X-ray 
spectrum must be produced through single Compton scat- 
tering. This requires very energetic comptonising leptons. In 
order to fit the RXTE data, the model requires an electron 
temperature ~ 3 MeV. 

We note however that such a combination of small size, 
very low optical depth and large temperature is physically 
impossible. Indeed, as discussed below in section 1431 in Cyg 
X-l, the large luminosity and small emitting region make 
the compactness larger than unity. Therefore achieving a 
very small optical depth and a temperature kT > m c c 2 is 
impossible due to the effects of pair production. Svensson 
et al. (1984) studied the equilibrium properties of a ther- 
mal pair plasma. His Fig. 6 shows that for a compactness 
of order of a few (like in Cyg X-l), the optical depth must 
be at least 10 times larger than in this jet model. We also 
made some simulations using the code of Belmont, Malzac 
& Marcowith (2008). Setting the optical depth of ionization 
electrons r p — 0.01 and a bolometric luminosity of 2xl0 37 
erg s~ emitted in a region of 30 Rg, we found a total equi- 
librium optical depth tt — 0.6 and a lepton temperature 
fcT e ~ 300 keV. The exact values of the equilibrium temper- 
ature and optical depth depend on the strength of magnetic 
field that, following M05, we assumed close to equipartition 
with the radiating electrons. The resulting X-ray spectrum 
was clearly different from both the observed and the ones 
computed by M05 neglecting pair production. 

Then we explored extensively the parameter space and 
found that it is not possible to have simultaneously a low 
optical depth and a temperature as large as the one required 
by the M05 model unless the size of the emitting region ro 
is larger than ~ 10 4 Rg- On the other hand if the emitting 
region is very compact ro < 10 the pair optical depth can 
reach unity, in agreement with usual Comptonisation fits, 
however the equilibrium pair temperature is then too large 
(> 150 keV) to reproduce the sharp cut-off that is observed 
in the hard X-ray spectrum around 100 keV. 

The parameters of the model of M05 are therefore in- 
consistent with the constraints from pair equilibrium. Once 
the effects of pair production are taken into account, it seems 
impossible to fit the high energy spectrum of Cyg X-l with 
this model. Finally we note that X-ray dominated jet models 
also appear to be ruled out by the comparison of the prop- 
erties of black holes and neutron stars (Maccarone 2005) . 

4-2.2 Hot flow models 

Our results are consistent with the popular idea that the 
X-ray emission is produced in some sort of hot accretion 
flow (e.g. Shapiro Lightman & Eardley 1976; Narayan & Yi 
1996). As we have shown, if the hot flow constituted the base 
of the jet (i.e. was being ejected to infinity), its ejection ve- 
locity would have to be very slow. Then, the ejection time 
could be comparable to, or even longer than the viscous time 
in the inner part of the accretion flow. If so the corona would 
be accreted before being ejected. It is therefore unlikely that 
a large fraction of the accreting material is ejected. Also for 
the reasons discussed in section 14.1.2 1 it is much more plau- 
sible that only a small fraction of the hot flow is loaded into 
the jet and ejected to infinity. In this framework, our con- 
siderations would be consistent with the coupled ADAF/jet 
model of Yuan, Cui & Narayan (2005). We stress however 
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that the accretion efficiency changes at most by a factor of 
~ 2.7 during the state transition (see section [2~Tj) and there- 
fore Cyg X-l cannot be strongly advection dominated in the 
hard state. In this context, radiatively efficient hot flow solu- 
tions such as the luminous hot accretion flow model of Yuan 
(2001) would be favoured. We also stress that in the X-ray 
emitting region of Cyg X-l, the proton temperature seems 
much lower than the predictions of standard ADAF models 
(Malzac & Belmont 2009). 

4-2.3 Dynamic accretion disc corona model 

Beloborodov (1999) and then Malzac, Beloborodov & Pouta- 
nen (2001) showed that the hard X-ray emission of Cyg X-l 
can be understood in terms of an accretion disc corona out- 
flowing with a midly relativistic velocity (/3 ~ 0.3) atop a 
geometrically thin, optically thick accretion disc. Since this 
velocity is comparable to the infered terminal radio jet ve- 
locity and also comparable to the initial launching velocity 
of some popular jet models, it is tempting to associate this 
dynamic corona with the base of the jet. However, the typi- 
cal parameters of a dynamic corona (f3 ~ 0.3) lie at the edge 
of the allowed parameter regime (see Figs [3] and [4} . Besides 
the extreme accretion efficiencies required, this would im- 
ply a very small jet section (<4 r g ) and may require the 
jet to decelerate so that /3oo < /3q. A far more likely situa- 
tion would be that the whole corona is indeed accelerated 
to relativistic speed (as in the dynamic corona model) but 
only a small fraction of its material actually escapes to in- 
finity. This could be related to the velocity of the gas being 
lower than the escape velocity close to the black hole (see 
e.g. Ghisellini, Haardt & Matt 2004). Indeed, within 20 R g , 
the escape velocity (/3 CSC = \Pzfr) is larger than 0.3c and 
most of the ejected gas would be accreted again. Whatever 
the origin of the X-ray emission (hot accretion flow or ac- 
cretion disc corona) we favour a situation in which only a 
small fraction of the X-ray emitting material is loaded onto 
the jet and ejected to infinity. 

4.3 Caveats 

Let us now consider the possible limitations to our results. 
First, our estimates of the jet terminal velocity and in par- 
ticular the lower limit obtained for this velocity relies on our 
assumption that the optical nebula of Cyg X-l is powered 
by the kinetic power of the jet. We note that the possibility 
that the nebula is actually a background supernova rem- 
nant cannot be fully ruled out (Russell et al. 2007). Even 
if the nebula is powered by the jet, it is possible that the 
energy of the jet is not dominated by kinetic energy (as in 
the case of a Poynting flux dominated jet for instance). If 
so, our analysis remains valid provided the j parameter is 
reduced accordingly. Then, the estimated jet velocity would 
be lower. 

Another possible limitation is that our estimate of /3oo 
relies on the absence of jet in the soft state. In fact, we can- 
not exclude that the HSS may produce a jet of very high 
Lorentz factor. Such a jet could transport away a large frac- 
tion of the energy and yet be unobservable because our line 
of sight falls outside the beaming cone (Maccarone 2005) or 
because of the absence of a shocked deceleration region. In 



our framework, this would be equivalent to reducing the ra- 
diative efficiency in the HSS and we would then infer lower 
jet velocities. 

However the good agreement with the estimates of the 
jet velocity obtained from the radio observations indicates 
that our assumptions are at least roughly valid. This sug- 
gests that the jet is not strongly dominated by Poynting flux 
and that any jet formed in the HSS does not take away a 
substantial fraction of the energy. We also note that even if 
this was the case, our constraints on the initial jet velocity 
would remain valid provided the radio estimates of the jet 
bulk velocity are correct. 

Because the electron to proton mass ratio is so low, we 
have neglected the contribution of leptons to the jet kinetic 
power. Taking into account the leptons would amount to 
divide our estimates of — 1 by a factor l+£m c /m p , where 
£ is the number of leptons per proton in the jet. We see that 
our results would be affected only if the jet composition is 
strongly dominated by electrons-positron pairs. Again, the 
agreement with the radio estimates for our typical set of 
parameters suggests this is not the case. 

We note however that an important contribution of 
electron-positron pairs to the optical depth r at the base 
of the jet would significantly affect the product joPoTo 
(see equation @) which would be increased by a factor 
£/(l+£m e /m p ). If pairs are important at the base, the initial 
velocity may be large. For a typical size of the X-ray emission 
of ro = 30R g and Lh — 2 x 10 37 the compactness parame- 
ter is larger than unity (I = 8) and pair production may be 
important. In order to investigate this possibility, we per- 
formed numerical simulations using the code eqpair (Coppi 
1992). For the parameters producing a hard X ray spectrum 
similar to that of Cyg X-l, we found that 1 ^ £ ^ 2, as long 
as the size of the X-ray emitting region ro > 10 (i.e. I < 24) 
and t > 1. Therefore if ro > 10 the effects of pairs would be 
weak, increasing the initial velocity by less than a factor of 
2. However if the X-ray emitting region is smaller, electron- 
positron pairs may dominate. If so the initial velocity could 
indeed be relativistic. 

Finally, in our calculations, we have assumed that the 
source of jet power was in the accretion disc. We have ig- 
nored the possibility that the jet may tap a significant frac- 
tion of its energy from the black hole spin (Blandford & 
Znajek 1977). If this additional source of energy is impor- 
tant, the net effect is to increase the apparent efficiency r/h 
in the LHS. Our results would be affected only if this pro- 
cess can lead to n^ > 1. We believe this is unlikely, specially 
if the black hole spin in Cyg X-l is as small as inferred by 
Miller et al. (2009) . 



5 CONCLUSION 

The jet power estimates of Russell et al. (2007) and the fact 
that the luminosity of the X-ray source does not increase 
dramatically during the state transition put some interest- 
ing constraints on the energetics of the accretion flow in 
Cyg X-l. Notably, the accretion efficiency cannot increase 
dramatically across the state transition from LHS to HSS. 
The jet bulk velocity must be relativistic (/?oo > 0.1) and, 
depending on the accretion efficiency, it is likely to be in 
the range 0.3-0.8. Then if the Thomson optical depth at 
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the base of the jet is larger than unity the initial jet veloc- 
ity must be very low. This is in contradiction with several 
jet models in which the X-ray emission is produced in the 
jet or its base. We also pointed out that the specific jet 
model of M05 appears inconsistent with the observed X-ray 
spectrum of Cyg X-l once the effects of pair production are 
taken into account. Finally, both hot accretion flow and out- 
flowing accretion disc corona models remain consistent with 
our results provided that only a small fraction of the X-ray 
emitting material is loaded on the jet. The X-ray emitting 
region and the jet therefore appear to be distinct although 
physically connected. 
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